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Abstract 

In this note we investigate U(N) gauge theories with matter in the fundamental and 
adjoint representations of the gauge group, interacting via generalized Yukawa terms of 
the form Tr [Q§ n Q]. We find agreement between the matrix model and the gauge theory 
descriptions of these theories. The analysis leads to a partial description of the Higgs 
branch of the gauge theory. We argue that the transition between phases with different 
unbroken flavor symmetry groups is related to the appearance of cuts in the matrix model 
computation. 
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1 Introduction 



Recently, Dijkgraaf and Vafa have proposed [Q, EL |3|] a method for computing pertur- 
batively the effective glueball superpotential of M = 1 theories with fields transforming 
in the adjoint and bifundamental representations of the gauge group. According to this 
proposal, the planar free energy of the matrix model whose potential is the tree-level 
superpotential of the J\f — 1 theory yields the effective superpotential of this theory. 

When fields transforming in the fundamental representation of the gauge group (quarks) 
are present, one only needs to include the planar free energy coming from diagrams with 
one quark boundary ||, ||. More explicitly, the gauge theory effective superpotential is 
proposed to be 



W cS (S, A) = N C S(1 - In — ) + N c ^l + N f F x=1 . (1) 



This prescription was successfully used to compare matrix model predictions with 
known gauge theory results for theories with massive and massless flavors, with M = 1 
and M = 2 supersymmetry @]-[10|. 



For theories with fields transforming in the adjoint representation of the gauge group, 
proofs that planar graphs are the only ones which contribute to the matrix effective super- 
potential were presented in Jll] (based on the analysis of superspace Feynman diagrams) 



and |12] (based on holomorphy and symmetries). The latter arguments were extended in 
T3| to the case of theories with fields transforming in the fundamental representation of 



the gauge group and it was shown that only planar diagrams with one (appropriately gen- 
eralized) quark boundary contribute to the gauge theory effective superpotential. Other 



interesting related work has appeared in |14 



The correspondence between gauge theories and matrix models has been pushed very 
far for superpotentials depending only on the adjoint fields. However, these checks have 
only been performed in the simplest cases of theories with fields transforming in the funda- 
mental representation of the gauge group. While it seems reasonable that the arguments 
of |T2| and |T3[ generalize to generic superpotentials, it would be interesting to perform 



some explicit checks, along the lines of [|J and []5J 

In this paper we work out the details of the matrix model and the gauge theory for the 
Q§ n Q coupling. With these results as a starting point, we then outline how a polynomial 
of generalized Yukawa couplings QP(§)Q can be analyzed. We find complete agreement 
between the matrix model with one boundary and the gauge theory. The Coulomb branch 
of such theories was discussed in detail in [[32|1 . However, the Higgs branch seems largely 
unexplored. The matrix model computations suggest a simple way for analyzing it. 

In the next section we use the matrix model to compute the effective superpotential 
of this theory. Since the adjoint field does not interact with itself, this superpotential is 
just the sum of a Veneziano-Yankielowicz piece (coming from the dynamics of the gauge 
field) and the free energy given by diagrams with one quark boundary. 

We find that the sum of these diagrams gives a free energy identical to that of a theory 
containing an adjoint and n quarks with regular Yukawa couplings giQi&Qi, where the 
coupling constants are proportional to the n roots of the unity. We then show that one 
can go from the second theory to the first by simply integrating out certain combinations 
of the n quarks until only one quark and the adjoint are left. 
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We then discuss the gauge theory origin of the matrix model results. The nonpertur- 
bative contribution to the effective superpotential of the theory with a QQ n Q coupling is 
hard to obtain by symmetry and holomorphy arguments. One might hope that integrating 
$ out might make things better, since only quarks will be left and the nonperturbative 
contribution to the superpotential would be of Afneck-Dine-Seiberg type |3B|]. Neverthe- 
less, after integrating out <3>, one is left with a "tree level" superpotential which contains 
the coupling constant to a negative power. Since this term does not have a well defined 
limit as g — > 0, one can no longer argue that this term cannot mix with the Afneck-Dine- 
Seiberg contribution. Therefore, one expects nonperturbatively generated terms which 
contain combinations of A and g. These nonperturbative terms cannot be easily found 
using analyticity and charge conservation. 

It appears therefore that in order to find the nonperturbatively generated contribution 
to the superpotential one has to find a theory (related to the theory of interest by inte- 
grating in and integrating out) where the nonperturbative contribution has a simple form. 
Fortunately, as the matrix result hints also, the theory with Yukawa coupling Tr [Q<$> n Q] 
can be obtained by integrating out n — 1 nontrivial linear combinations of quarks in a 
theory with an adjoint, n quarks, and Yukawa couplings Z)™ =1 e 2ml ^ n Tr [Qi&Qi\- At this 
stage one can integrate out the adjoint field and find a theory with n quarks and interac- 
tions of the form Tr [(J2i 9iQiQi) 2 ]- F° r this theory one can use the usual holomorphy and 
charge conservation arguments [[34] to show that the only nonperturbatively generated 
superpotential is the Afneck-Dine-Seiberg one. 

Once we have the full gauge theory effective potential it is not hard to integrate out all 
the fields and relate the resulting effective superpotential W e g(A, rrii, gi) with the matrix 
model computation. The Higgs branch of the original theory can also be analyzed. 

We emphasize that his effective superpotential is the same, regardless of the order in 
which one integrates out the fields^], and thus regardless of the hierarchy of the m r This 
is quite obvious from the matrix model perspective: by summing all (planar) Feynman 
diagrams one obtains the same function of mass parameters, regardless of their hierarchy. 

One can also see this in the gauge theory. The physical mass of a field depends both 
on its superpotential mass parameter m, and on the Kahler potential. Since by adjusting 
the latter any hierarchy can be achieved (regardless of the magnitudes of the m's), and 
since this adjustment does not affect the superpotential, it follows that the final result is 
independent on the mass hierarchy and on the order one integrates fields out. 

Therefore, the effective superpotential W eS (A, m^, gi) one finds after integrating out 
all the fields is the same if the theories we start from can be related to each other by 
integrating in or integrating out. In the case of the theory we discuss, this effective 
superpotential is most easily obtained by considering the related theory with n quarks 
and simple Yukawa interactions, finding its nonperturbatively generated superpotential, 
and integrating out all the quarks. This computation appears in the last section of this 
note. 

Note Added: When this work was near completion we received the preprint [|Kj 
which, while having a different focus, overlaps with the technical details of our work. 

4 We are grateful to Eric D'Hoker for pointing this out to us. 
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2 Matrix model 



As we recalled in the Introduction, the part of the matrix model free energy J-'(S) which is 
related to gauge theory via extensions of the DV prescription is computed by summing all 
the planar Feynman diagrams with one quark boundary. If the quark-adjoint interaction 
is of the form Q^Qi, the necessary combinatorics of the Feynman diagrams was described 
in |4], [35]]: a diagram with 2k vertices comes with a factor of T^ky from the exponential; 



then, the different ways of contracting the quarks produce a factor of (2k — 1)!; the number 
of ways of connecting, in a planar manner, the boundary points with adjoint propagators 
give rise to a factor , fc ^jj fc! . Since each diagram contains 2k quark propagators, k adjoint 

propagators, and 2k vertices, it is multiplied by (^ M ^ m -2 ^ • Finally, the external flavor 

loop gives a factor of Nf, while the k + 1 color loops give a factor of S k+1 , where S is the 
't Hooft coupling of the matrix model, and becomes identified under the correspondence 
with the gauge theory glueball superfield. 

Thus, the free energy contribution of these diagrams is f| 

^ = ~ Nf § frw sk+1 (2) 

In the case of diagrams with a g n Qi& n Qi interaction, it is necessary to distinguish 
between the case of odd and even n. If n is odd, only diagrams with even numbers of 
insertions contribute. The combinatorial factors coming from the quarks are unchanged. 
Nevertheless since now there are nk $ lines, there will be ( rafc +")^ fc )i different ways of 
connecting them, and the overall power of S will be nk + 1. Thus 

7r°dd _ at (2n£; - 1)! qnk+ i ( g\ \ , . 

If n = 2p is even, diagrams with any number of insertions contribute. The free energy 
is ^ 

= -2vN f T (2pfc- 1)! SP k+i ( 9n \ (4) 

One might have also expected extra factors of n\ coming from the different orderings 
of the $'s originating from one vertex. However, since the $'s are matrices, one cannot 
interchange them at a vertex because this would make the diagrams nonplanar. Thus, 
the diagrams give the same answers as if the n $'s originating at one interaction vertex 
were separated by tiny propagators of some auxiliary quarks. This is depicted in Figure 
1, and is a hint toward the equivalence of our theory to a theory with n quarks and simple 
Yukawa couplings, equivalence which will be discussed in the next section. 

If one introduces the new variable p, such that pn = n where 

k = 2 if n is even, 

k — 1 if n is odd , (5) 

it is not hard to see that the free energy for both even or odd n is given by: 

T x =\ = -KpNfKp, (6) 
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Figure 1: Equivalence between Yukawa and generalized Yukawa couplings. 



where K p is the value of the sum in 



The radius of convergence of the sum JC P can be easily found to be S c = 4-. For S < S c 



the sum is: 



lP -i 



1 

2 + 4( 



p<yS 



1 -A(-fvaS- 1 



lni(l + yi-4(-) 2 ^) 



(7) 



and a r 



9n 



For S > S c the sum is divergent, and one must find its value by analytical continuation 
[fTOfl . Since (0) contains square roots, one expects branch cuts in the complex S plane 
starting from the points where the square roots become zero and ending at infinity. 

Moreover, since we have n square roots, each comes with a choice of branch. Therefore, 
for 5* outside the radius of convergence, the sum (0) has 2 n values, depending on the choice 
of branch for each square root. Furthermore, since the sum ([?]) applies separately for each 
of the flavors, there will be a total of 2 nN f branches for the free energy. The result is: 



/Cp(e) 



1 



l 

2 + ^ 



'paS . 



4( 



In 



1 + ei)/ Wl-4 



~>aS 



with eij = ±1 with / = 1, . . . , n and / = 1, . . . , Nf . 

For the example discussed in |l(J, the choice of branch in the matrix integral was 
matched in gauge theory with the choice of roots of a certain the second order equations. 
However, one can choose the parameters of the theory such that all relevant values of 
S lie inside the radius of convergence. We limit ourselves to showing agreement in this 
regime. We will show that the convergence of (0) for S < S c has a precise meaning in 
gauge theory. The branch structure and the corresponding phase structure of the theory 
|IIJ follows from the careful analytical continuation of our results. 

To evaluate the superpotential ([!]) at its critical point we begin by constructing a single 
formula which covers both cases k = 1 and k = 2. Since the set of odd roots of unity goes 
to itself when squared, it is not hard to see that the free energy can be expressed solely 
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in terms of n: 



n—1 

1=0 



Ae 2m ~aS 



1 - Ae 2m ^aS - 1 



ln^(l + y/l-4e 2 ™%aS) 



(9) 

Since all roots of order n of unity appear in the above expression, it is clear that all phases 
in the definition of a below equation (|7|) are equivalent. 

It is quite easy to see that each term in the sum in equation (|9|) reproduces the 1- 
boundary free energy of a theory with a single quark and a regular Yukawa coupling. 
The ratios of the coupling constants of these theories are n-th roots of unity. In the next 
section we will show how this comes about in the gauge theory. 

Using (^), the critical points of the effective superpotential ([]]) are given by the solution 
of the equation: 



In 



S 



n-l 



1 



N f A 3 



n-l 



£>-(! + yi-4:e 2m taS) = In J] -(1 + ^l-Ae 2 ^aS) 



(10) 



7 ^ i=o 
which can be trivially transformed into: 



1=0 



N c " 1 / 

y -f =]T -(l + V / l-4e 2 ^/?y) , 
i=i A 

with y = T5 and (3 = aA 3 . 

Then, the values of the superpotential at its critical points are given by 

1 



w\ 



7? 

N c --N f + N f J2 

A i_i 1 I , H A ^2tti- 



i=i 1 + VI -4:e 2m ~py. 



y 



(12) 



where we have to replace y by a solution of ([TT|) . This is the result that we will compare 
with the gauge theory predictions. We stress that this equation is valid only in the limit 
of small /3, i.e. far from branch points of the series (0). 

For generic (3 the critical values of the superpotential are given by 



W\ 



A 3 



N 



f n 



f=l 1=1 1 



l-4e 2m -{3y 



(13) 



where y is given by an equation similar to (JTTj) , except that the square roots are dressed 
with e coefficients. 



3 Gauge theory 

Let us begin the gauge theory analysis by showing that the generalized Yukawa coupling 
Tr [Q<& n Q] is equivalent to a set of n ordinary Tr terms whose strenghs differ by 

n-th roots of unity. To this end we notice that, by starting with the superpotential 

W tree = -M Tr $ 2 + m £ Tr [Q^] + g ^ Tr [Q^Q l+1 ] with Q n+1 = Q x (14) 
1 i=i i=i 
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and integrating out Qi and Qi for all i — 2, ... n we recover 

W trec = X -M Tr $ 2 + mTr [QQ] + gTt [Q$ n Q] 

provided that m™ -1 ^ = g. 

Equation flU]) can be rewritten as 



(15) 



-MTr$ 2 + m^Tr[g,4] + 



(n-l)/n 



(16) 



with the matrix V being given by 



V 








. 


. 


1\ 


1 





. 


. 








1 


. 


. 





Vo 





. 


. 1 


0/ 



(17) 



Since we chose the mass matrix of the nNf quarks to be proportional to the identity 
matrix, it is clear that the Yukawa and the mass terms can be simultaneously diagonalized. 
Noticing that eigenvalues of V are given by the roots of unity, we can immediately rewrite 



(0) as 

* ' tie 



-MTr $ 2 + Tr +X>tf lWQj$Qi] with = e 



2m ^ 



1=1 



We have thus shown that the coupling Tr [Q§ n Q] is equivalent to a set of diagonal Yukawa 
couplings whose strengths differ by roots of unity 

Unlike equation (|14|), the superpotential (|18D is invariant under global SU (Nf)® n trans- 
formations. In constructing the Affleck-Dine-Seiberg superpotential it is useful to think 
of (PD particular case of 



^ treo = -MTr$ 2 



mTr [QQ] + Tr [GQ&Q] 



(19) 



where now Q are in SU(nNf) and G is a diagonal matrix. The SU(nNf) invariance if 
broken either by having different coupling constants or by having different masses for the 
n sets of Nf quarks. 

Since we are interested in comparing this gauge theory to the matrix model of the 
previous section, we first integrate out the adjoint fields. The effective superpotential is 
given by: 

W tICC = mTr [QQ] - 7^ Tr [GQQGQQ] . (20) 

We will later identify the coupling constant matrix with the one given by equation fll8D , 
but we will derive the formulae for an arbitrary diagonal G = diag(gi, . . . , g n ). 

To this tree level superpotential we have to add the nonperturbative contributions. It 
is not hard to see that they are given by the ADS superpotential for Q and Q. Indeed, 
in the limit of vanishing m and G this is the only possible term. Demanding analyticity 



G 



as well as preservation of the symmetries leads to the conclusion that no corrections are 
possible. 

Next we want to integrate out all quarks. The easiest way to find the result is to rewrite 
the above superpotential in terms of mesons and notice that, since the mass matrix as 
well as G are diagonal, all off-diagonal components of the mesons are constrained to 
vanish. Writing the remaining components of the meson field as Xa = QiQi = a^l/Vy the 
remaining effective superpotential is: 



W eS = N f mJ2xi-N f J2 a i x l + ( N c ~ nN f ) 



i=i 



i=l 



\3N c -nN f 

rij=i %i f 



(21) 



where a% = 4h- Minimizing this superpotential gives 



3N c -nN f n 



mxi - 2a;x- - A Nc ~ nN f JJ 







(22) 



i=i 



To compare with the matrix model, it is useful to make the following change of variables: 

N f 



nN c -nNf 



(23) 



in terms of which the equations of motion can be rewritten as: 



y% - PiVi - n Vj 



N c — n,Nj- 



= 



(24) 



j'=i 



where we defined 



m 



jy N c -nNf 



m 



(25) 



and we recall that % — 1, . . . , n labels the different types of quarks. This gives a system of 
n coupled non-linear equations. To proceed it is helpful to introduce a "radial" variable: 



3=1 



(26) 



The equation of motion becomes 

Vi-PiVi-V = Q ■ (27) 
We can solve for each individual yi in terms of the couplings and the radial variable: 

1 



Vi = 



2A 



■1 + £^1-4/3, 



where e« = ±1. Each choice of e-s one finds an equation for y: 



7 =nzr 

i=i 1 



-2A 



+ e iy /l - iPiy 



ri2A 



1 + gjVj - A(3 iy 

my 



(28) 



(29) 
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Bringing factors of y together, we get the following equation: 



n 1 / 

JI-(l + e 4> /i_4Ay) 
1=1 z 



(30) 



This algebraic equation can be solved numerically for various values of n and the couplings. 
Once one has a solution of this equation, one can obtain the y^s from equation fl28|) for 
each of the 2 n choices of 6j. 

We next compute the effective superpotential at the minimum: 



W\ 



crit 



lN f mjTxi + \(2N c -nN } ) 

i=l ~ 



^3N c -nN f 



Nc—nNf 



nNj 3N c -nNj, 



i n i n 

^Nf^Vi + ^Nc-nN^Hy, 

i=l " i=l 



(31) 



In fact, the superpotential can be written in terms of the variable y alone by using equation 



W\ 



A 3 



A 3 



1 1 n -1 

-{2N c -nN f )y+-NfY { — ( -1 -r,- v <l- 1 ?,// 



N c - 



-Nf + NfY* ; 



(32) 



where the scale L is defined by 



A 3 



"A" 



(33) 



which is the correct relation between scales when all nNf quarks are integrated out. 

So far the discussion was for general Pi. Now, when all of the couplings /?/,/ = 1, n, 
are different, the number of solutions to the system ( P?[ ) is 2 n times the number of solutions 
of equation (P?p. If we chose a more general form of the meson, Xu = diag(xj, . . . , xf f ), 
the number of different vacua would increase to 2 nN f times the number of solutions of the 
apropriately modified equation (^). This matches the total number of branch cuts in the 
matrix model computation. Furthermore, it is easy to see that specializing the coupling 
constants g\ to the ones implied by the equation fli~8|) we find 



e 2 ^ a 



m 2 l n 2M 2 



Pi 



~(3 



(34) 



Therefore, the gauge theory and the matrix model results match in the region where 
the series (0) develops branch cuts. However, since for y(3 < | the series leading to (|) is 
convergent, all e^/ coefficients are fixed to unity, while in gauge theory the choice of signs 
€ij seems to persist for all values of y. 

To understand the solution of this apparent discrepancy we should note that the ap- 
pearance of different branch points leads to a spontaneous breaking of the U(Nf) flavor 
symmetry in a theory with Nf identical quarks. In the regions of parameters where both 
€{j = +1 and €ij = — 1 are allowed, there exist vacua with broken flavor symmetry. 
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However, one expects that, as the coupling constant is reduced, the flavor symmetry will 
be restored. Indeed, by taking the small fa limit on equation ( GBp we find that, for any i, 
the solution corresponding to 6j = — 1 moves off to infinity while the one corresponding 
to 6ij = +1 remains at finite distance. Thus, in the small coupling limit, only the choice 
6i = +1 is allowed. 

This has definite meaning in the matrix model computation. In the small coupling 
limit the radius of convergence of the series (]?]) becomes very large. Therefore, all branch 
points move off to infinity and there remains a unique choice for the free energy T x= \. 

We are therefore led to interpret the appearance of branch cuts in the series (0) as the 
matrix model version of transitions between domains of the Higgs branch with different 
flavor symmetry. 



4 More generic quark-adjoint interactions 

We can generalize the techniques we developed for Tr [Q& k Q] interaction to study theories 
with superpotentials of the form Tr [QP(§)Q], where -P($) is an arbitrary polynomial of 
finite degree. 

The matrix prescription is obvious. One must compute quark contribution to the 
free energy by summing all the one boundary Feynman diagrams with arbitrary types of 
insertions. The vertices are given by the monomials appearing in -P($), and the quark 
and adjoint propagators are the inverse of their masses. 

While a formal expression for the free energy of the matrix model can be written for a 
generic -P($), it does not seem of any particular use. Let as only illustrate the procedure 
by considering a theory with an interaction term of the form 

Tr [Q( 9l ^ + g 2 ^ + gs$ 2p3 )Q] , (35) 

where we have chosen the powers of $ to be even in order to keep the counting easy. 
A general diagram containing hi vertices of type i comes with a factor of 1/(A; 1 ! ^2^3!) 
from the expansion of the exponential; the different ways of contracting the quarks give 
a factor of (ki + k 2 + k 3 — 1)!; the different ways of connecting the boundary points with 



non- intersecting adjoint propagators give [35] a factor 



(2k lPl + 2k 2P2 + 2k 3 p 3 )\ 

(36) 



(km + k 2 p 2 + k 3 p 3 + + k 2 p 2 + k 3 p 3 )\ 



Since each diagram contains k\ + k 2 + k 3 quark propagators, km + k 2 p 2 + k 3 p 3 adjoint 
propagators, and k±, k 2l and respectively k 3 vertices, it is multiplied by 

fci &2 &3 

<Jl '.)> U:\ 



T\ / fkipi+k2P2+k3P3 rr ,k 1 +k 2 +k 3 
1V1 <S> " l Q 



Finally, the external flavor loop gives a factor of Nf, while the (k\p\ + k 2 p 2 + k 3 p 3 + 1) 
color loops give the appropriate power of the glueball superfield. Thus, the 1-boundary 
contribution to the free energy of this model is simply 

kl ,k 2 ,H=o {Ki\k 2 \k 3 \) 
k l +k 2 +k 3 ^0 
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x (2fc 1 p 1 + 2fc 2 p 2 + 2fc 3 p 3 )! (S_ 

(kipi + k 2 p2 + hP3 + + hP2 + hp 3 ) ] - yMtj, 

For odd monomials the above sums have to be only over combinations of terms which 
give an even number of $ propagators. It is quite easy to see that in general such sums 
are hard to compute explicitly. 

Nevertheless, it is possible to obtain this free energy by using the vertex splitting 
procedure we used in section 3. A vertex of the form g n Ti [Q<& n Q] can be thought of as 
arising from a theory with n quarks and off diagonal Yukawa interactions gTr [Q n §Q n V n } 
of the type (|TBD by integrating out the auxiliary quarks Qi. 

Therefore, the matrix model of the theory with a polynomial interaction Tr [QP(§)Q] 
can be related to that of a theory with interactions linear in $ if one introduces n — 1 
auxiliary quarks for each monomial g n Tr [Q$ n Q] coming from P($). One obtains an off- 
diagonal interaction matrix whose dimension is the sum of the powers of the monomials 
in P($) minus the number of monomials plus one. Since the auxiliary quarks have 
the same mass as Q 0, one can diagonalize the interaction matrix and obtain a theory 
with diagonal Yukawa interactions AjTr [QjQQj], where the Xj are the eigenvalues of the 
interaction matrix. 

As in the case discussed in section 2, the partition function will be the sum of 1- 
quark regular Yukawa partition functions with couplings \j. The only difference from the 
equation (|7|) will be that the couplings will not be proportional to the roots of unity, but 
will have a more complicated form. 

To treat this theory correctly in the gauge theory, one must perform the same steps, 
by integrating in the auxiliary quarks and obtaining a theory with only linear couplings 
of the adjoint field $. One can then integrate out $ and obtains a theory with only 
quarks. For this theory it is possible to determine that only the regular Afneck-Dine- 
Seiberg superpotential is generated nonperturbatively; one can then integrate out all the 
quarks and obtain a superpotential which can be related to the matrix one. 

It appears therefore that since integrating in and out work identically on the two 
sides, the equivalence of the matrix result and the gauge theory result is ensured by the 
equivalence of these results for theories with quarks with equal mass and different Yukawa 
couplings. This equivalence is obvious from the computations in sections 2 and 3, and 
follows also from the results of |TU| by making particular choices for the mass parameters 
and rescaling the quarks. 

5 Conclusions 

We have investigated an U(N) gauge theory with adjoint and fundamental matter in- 
teracting via a coupling of the form Tr [Q<& n Q]. We have solved the matrix model and 
found the exact low energy effective superpotential. This effective superpotential is iden- 
tical to that of a theory with n quarks minimally coupled to $, with coupling constants 
proportional to the n'th roots of unity. As expected, these two theories are related by 
integrating in/out n — 1 quarks. 

On the gauge theory side we argued that in order to determine unambiguously the 
nonperturbatively generated contribution to the superpotential one needs to first integrate 

5 As we explained, this does not affect integrating them out. 



fel+fc2+fc3 

m o 
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in these auxiliary n— 1 quarks, obtain a theory with minimal couplings between the quarks 
and the adjoint field and then integrate out $ to obtain a theory with n massive 
quarks and a quartic tree-level superpotential. One can then use standard holomorphy 
and symmetry arguments to argue that the only nonperturbative superpotential in this 
theory is the Affleck-Dine-Seiberg one. By integrating out all the quarks we obtained the 
low energy effective superpotential of our theory, and we found it agrees with the one 
computed in the matrix model. 

We also described a method to investigate theories with more complicated adjoint- 
quark couplings, of the form Tr [QP(<&)Q], where P($) is a generic polynomial of finite 
degree. We illustrated this technique by writing down the matrix free energy for a poly- 
nomial P built out of three monomials of arbitrary even powers. While the generalization 
is straightforward, the perturbation theory is difficult to resume. It is also possible to 
further generalize this discussion by adding an arbitrary superpotential depending only 
on the the adjoint field. 

We then presented a method to solve these theories by relating them to gauge theo- 
ries with many quarks but only minimal couplings A«Tr [Qi<&Qi\. When the polynomial 
contains just one monomial of order n (this is the case discussed in the first two sections 
of this note), the \ are proportional to the n'th roots of unity. For more complicated 
polynomials, the \ are the eigenvalues of the quark interaction matrix. 
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